Supersolidity and phase diagram of softcore bosons in a triangular lattice 



Jing-Yu Gan, 1 Yu-Chuan Wen, 2 and Yuc Yu 3 

1 Center for Advanced Study, Tsinghua University, Beijing, 100084, China 
^Institute of Theoretical Physics and Interdisciplinary Center of Theoretical Studies, 
Chinese Academy of Sciences, Beijing 100080, China 
s Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100080, China 

We study the softcore extended Bose Hubbard model in a two-dimensional triangular lattice by 
using the quantum Monte Carlo methods. The ground state phase diagram of the system exhibits 
a very fruitful structure. Except the Mott insulating state, four kinds of solid states with respect 
to the commensurate filling factors p — 1/3, 2/3 and p = 1 are identified. Two of them (CDW If 
and CDW III) are newly predicted. In incommensurate fillings, superfluid, spuersolid as well as 
phase separation states are detected . As in the case for the hardcore bosons, a supersolid phase 
exists inl/3<p<2/3 while it is unstable towards the phase separation in p < 1/3. However, this 
instability is refrained in2/3<p<l due to the softening of the bosons and then a supersolid phase 
survives. 
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Introduction A quantum matter state with non-zero or- 
der parameters characterizing both solid and superfluid 
phase, so-called supersolid (SS), was proposed by Pen- 
rose and Onsager fifty years agopj. In 1970s, it was 
speculated that the large zero-point quantum fluctuation 
in solid He may induce an off-diagonal-long-range-order 
0. The recent observation of non-classical rotational 
inertia (NCRI) in porous Vycor or in bulk solid 4 He re- 
ported by Kim and Chan [3 have greatly revived the in- 
terest in studying this new matter state although it is 
still controversial whether the mechanism suggested by 
Refs. @ works in explaining the NCRI 0]. 

Comparing to the complication of the 4 He, the di- 
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FIG. 1: (Color online) The ground state phase diagram for 
the extended Bose Hubbard model in the two-dimensional 
triangular lattice(£7 = 30t). SF is standing for the superfluid 
phase; MI for the Mott insulator phase. The p = 1/3 solid 
means one of three sets of sublattices is occupied with one 
atom per site. In the inset, the three different solid states 
are illustrated. In CDW I at p = 2/3, two sets of sublattices 
are occupied with one atom per site. In CDW II at p = 2/3, 
one set of sublattices is occupied with two atoms per site; In 
CDW III at p — 1, a single set of sublattices is occupied with 
three atoms per site. 



lute cold atom gas loaded in an optical lattice |FJ is a 
more promising avenue to success the SS, especially when 
the long-range interacting atoms were already cooled 
to a Bosc-Einstcin condensate These systems may 
be well described by the Bose-Hubbard model Q and 
its extended version. The simple Bosc-Hubbard sys- 
tem docs not allow a SS phase for the lacking of the 
interaction-based long range correlation. It was also 
checked that there is no a SS phase in the hardcore 
Bose-Hubbard model with a nearest neighbor repulsion in 
two-dimensional square lattice due to the solid-supcrfluid 
phase separation(PS), while a strip supersolid presents 
if the next nearest neighbor interaction is switched on 
||. The lattice frustration may stabilize a supersolid 
against the PS. This has been shown in the hardcore 
Bose-Hubbard model in a triangular lattice around the 
half filling 0, Hfl fill whereas the supersolid order is 
not favored in more complicated frustrated lattices such 
as Kagome lattice The next nearest neighbor in- 

teraction in triangular lattice can also lead to a strip 
supersolid in triangular lattice 0| ■ 

Besides the next nearest neighbor interactions and lat- 
tice frustration, softening of the on-site interaction may 
also stabilize the supersolid phase because more than one 
atoms per site arc allowed, which has been confirmed 
either for a two-dimensional square lattice |l "J| or one- 
dimensional lattice 0]. 

Model and explanation to results As we have seen, the 
physics of the supersolid phase is determined by the in- 
terplay among the next nearest neighbor interaction, the 
frustration of lattice as well as the softening of the on- 
site interaction. In this paper, we would like to check the 
cooperating effect of latter two. The model we study is 
the extended Bose-Hubbard model in a two-dimensional 
triangular lattice, which reads 

H = -t^2(a]a j + a]ai) + — 1 ^2ni{n i -l) 

(ij) 1 



2 




(a) (b) <c) 



FIG. 2: (Color online) The p = 1/3 solid doped by holes, 
(a) Placing the doped holes into one line costs no additional 
potential energy; (b) shifting the lower part of the lattice 
introduces a domain wall without any potential energy cost; 
(c) the doped holes (atoms) can hop freely across the domain 
wall. 
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FIG. 3: (Color online) The average atom density as a func- 
tion of the chemical potential for V = It, lOt and 12i, 
respectively (U = 30i). The system size is L = 12. 

where <xj( a ») is the creation(annihilation) operator of 
bosonic atom at site i; rii = a\a,i is the occupation num- 
ber; p is the chemical potential, (ij) runs over nearest 
neighbors. U and V represent the on-site and nearest 
neighbor repulsive interactions, respectively. Our central 
result is presented in Fig^ the ground state phase dia- 
gram. The quantum Monte Carlo(QMC) simulation with 
the stochastic series expansion(SSE) 0,0] was applied. 
We use t as the energy unit and focus our calculations on 
U = 30*. 

This phase diagram can be qualitatively explained by 
strong-coupling arguments 0, 0|. When U — » oo, the 
Hamiltonian (JJJ reduces to the hardcore one. The phase 
diagram is symmetric about p = 1/2, i.e., p/V = 3, 
because of the particle-hole symmetry [ij. The CDW II 
and CDW III as well as the upmost SS are ruled out. The 
phase transition from the solid to the SF is the first order 
one for the solid-SF PS appears at the phase boundary. 
Finite U breaks the particle-hole symmetry and results 
in other two solid phases and one more SS phase. The 



phase diagram with t = is quite easily understandable. 
A meaningful chemical potential is p > —U/2 otherwise 
p < 0. Raising p, there would be several kinds of solid 
phases with a common commensurate wave vector Q = 
(47r/3, 0). The first one is the p = 1/3 solid in which one 
set of sublattices is occupied with single atom per site. At 
p = 2/3, there are two kinds of solids depending on the 
ratio V/U. For small values of nearest neighbor repulsion 
Viy/U < 1/3), two sets of sublattices are occupied with 
single atom per site (CDW I), while for V/U > 1/3, one 
set of sublattices is occupied with two atoms per site 
(CDW II). For p = 1, the solid phases involved are the 
MI and the other one. When V/U > 1/3, one set of 
sublattices is occupied with three atoms per site (CDW 
III). If V/U < 1/3, the MI, i.e., single atom per site in 
the whole lattice, is the most stable. For p > 1, these 
commensurate solid phases are repeated in filling factors 
with an integer added. 

Based on this strong coupling limit understanding, we 
can describe the phase diagram we obtained. Comparing 
to the hardcore boson phase diagram in ref . , the p = 
1/3 solid lobe change slightly with only the up-bound 
lowered a little due to quantum fluctuations. When holes 
are introduced into the p = 1/3 crystal, the domain- wall 
proliferation mechanism, i.e., the additional holes hop 
freely across the domain wall, raises the kinetic energy 
linearly in t (see Fig . |2J and then excludes the possibility 
of a SS phase This leads to a PS between p= 1/3 

solid and a uniform superfluid with p < 1/3 [lflj . 

When atoms are introduced into the p — 1/3 crys- 
tal, the domain-wall proliferation mechanism does not 
work because formation of a domain wall would cost ex- 
tra potential energy. Thus, a SS phase is there right 
above this solid lobe, similar to that in the hardcore bo- 
son case. However, it completely separates the p = 1/3 
and p = 2/3 solid lobes. Furthermore, the latter shrinks 
much and splits into CDW II and CDW I in the V' 1 
direction due to the ratio V/U is lowered. 

Doping atoms into the p = 2/3 solid, the system be- 
haviors differently in terms of the ratio between V and 
U. Following an analysis in Ref. ^5} for the forma- 
tion of a SS phase above the half-filling in a square 
lattice, we argue that a stable SS exists for p > 2/3. 
For V/U < 1/3, CDW I is doped. Since the lattice 
sites are either single occupied or empty, the energy cost 
for an additional atom fills in empty or single occupied 
site is Eq = QV — p or E[ = U + 3V — p, respec- 
tively The energy of a single atom delocalized between 
these two cases is E 1 = E T Q + A 1 - ^/A 72 + (6t) 2 with 
A 1 = (U- 3V)/2. Similarly, for V/U > 1/3, i.e., doping 
in CDW II, an additional atom fill either in an empty 
or double-occupied site. The corresponding energy cost 
is either Eq 1 = 2U — p or E% = 6V — p, respectively. 
The kinetic energy of a single particle delocalized between 
these two cases is E u - E^ 1 = A 11 - ^A //2 +(6i) 2 , with 
A 11 = (3V - U)/2. In both cases, when A ~ t, the 
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FIG. 4: (Color online) Static structure factor S(Q) and superfluid density (p s )as functions of atom density p for V = 7t, lOt 
and 12t, respectively (U = 30t). Q = (47r/3, 0) for triangular lattice. 



kinetic energy E — Eq is proportional to t. Therefore, 
the domain wall formation was blocked and the system 
is stabilized against the PS As a result, these doped 
atoms form a SF on the top of the p = 2/3 solid back- 
ground, which is thought as a SS phase. Notice that the 
above analysis is not applied to U — > oo because multiple 
occupation is not allowed. The domain wall forms if dop- 
ing atoms into the p = 2/3 solid and then a PS happens 

For p = 1, the ground state is a solid with one set of 
sublattices occupied with three atoms per site (CDW III). 
Doping holes into this state, these additional holes move 
under the solid background with the effective hopping 
t* ~ t 2 /(7V — 2U). The kinetic energy gain remains 
quadratic in t and small. Therefore, the SS state with 
p > 2/3 is stable. However, our numerical result does 
not fully confirm this argument as we shall see later. 
Quantum Monte Carlo simulations We now portray our 
QMC simulations. In order to characterize different 
phases, we study the static structure factor S(Q) with 
Q = (47r/3, 0) and superfluid density p s , 

S(Q) = ^e i ^-^{n i n j ), p s = ^ (2) 

ij 

where W is the winding number fluctuation of the 
bosonic would lines [2(|; f3 = 2L is the inverse of fic- 
titious temperature and N = L x L is the lattice size. 
In our calculations, we use L = 9, 12 and 18. A solid 
phase is characterized by a non-zero S(Q) and p s = 
and a SF phase by S(Q) = and p s ^ 0; a SS phase is 
depicted by non- vanishing both of S(Q) and p s . Solids 
are categorized by their S{Q) values. 

Fig. |31 shows the atom density p varying as the chemi- 
cal potential p for V = It, lOt and 12i, respectively([/ = 
30t). It is seen that for each V, there is a jump from 



p < 1/3 to p = 1/3 which indicates a first order phase 
transition. In the grand canonical ensemble, a jump in p 
is a token of a PS region. For a small V, say V = 7t, only 
p = 1/3 solid plateau is observed, i.e., there is no p = 2/3 
plateau. This is different from that in the hardcore boson 
phase diagram for in that case there are both the p = 1/3 
and p = 2/3 solid lobes symmetrically. Obviously, this is 
because the particle-hole symmetry is broken by the fi- 
nite U. As V is raised, say V = lOt and 12t, the p = 2/3 
solid plateau appears and is gradually widened. Such an 
observation of the p = 2/3 solid plateau indicates the 
shrinking of the p = 2/3 solid lobe as shown in Fig. ^ 
In the regions 1/3 < p < 2/3 and 2/3 < p < 1, both 
of them include a SS phase. For V = 12t, another solid 
plateau at p = 1 emerges at larger p, while this plateau 
is not observed for V — 7t and lOt. 

In Fig. we show S(Q) and p s vary as functions of 
p. At p = 1/3, for all V we considered, the superfluid 
density p s — > while S(Q)/N ~ 1/9 which is an exact 
magnitude of the static structure factor of the p = 1/3 
solid state. For p > 1/3, the results are severely depen- 
dent on the magnitude of V. For a small, say V = 7t, as 
p increases, we observe, in turn, a p = 1/3 solid plateau, 
a SS phase in the region 1/3 < p < 0.4 and a SF after 
p > 0.4. No CDW III solid state is observed at p = 1. 
For an intermediate value of V{= lOi), both S(Q) and p s 
are finite and non-monotonous as p raises. A SS phase 
expands to the whole region from p = 1/3 to p = 2/3. 
At p = 2/3, p s -> while S{Q)/N - 0.1 which is the 
magnitude of the static structure factor of CDW I. With 
further increasing of p, there is another SS in between 
2/3 < p < 0.8. After p = 0.8, S(Q) -> 0, i.e., the sys- 
tem is in a pure SF phase. Again, no CDW III exists at 
p = 1. For large V, say V = 12t, p s — > at p = 2/3 while 
S(Q)/N - 0.45. This is CDW II solid state. As p in- 
creases further from p = 2/3, S(Q) raises monotonously 
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FIG. 5: (Color online) The ground state phase diagram for 
two-dimensional extended Bose Hubbard model in the trian- 
gular lattice in the V — p plane at U = 30t. 

and both p s and S{Q) in nonzero till p = 1. Thus, in 
2/3 < p < 1, the ground state is a SS phase. At p = 1, 
the ground state is another kind of solid state with p s — > 
and S(Q)/N ~ 1, i.e., CDW III. 

Phase diagrams and conclusions We scan V from 3i to 
12i with a wider step. This sketches the ground state 
phase diagram of the system . Fig. and Fig. 03 are 
its projections to the p—V^ 1 and V — p planes, respec- 
tively. The PS phase is explicitly shown in Fig. [S] and 
its width grows as V is strengthened. On the left hand 
sides of the CDW II and CDW III in Fig. we use dot 
lines labelling two possible PS regions. These PS states 
are neither supported by our early analysis nor confirmed 
by our numerical data. However, our data could not ex- 
clude these PS. What we see is although there seems no 
discontinuity in p of p, e.g., for V = 12t in Fig. [3] the 
change of p may be very sharp and for some specific mag- 
nitudes of V it is near a right-angle. We squint towards 
there are no these PS states but simulations with larger 
lattice size and lower fictitious temperature are required 
and confirmations with other numerical methods, e.g., a 
canonical ensemble calculation, are necessary. No such 
an ambiguity for the hardcore bosons because there are 
no CDW II and CDW III states 0j. 

At p — 1, for quite small values of V, say V = 3t 
and At, the MI phase( p s — > and S(Q) — > 0) exists. 
There is a SF region in between the MI and CDW III 
solid states and this region is suppressed as U becomes 
strong. Comparing with the hardcore boson phase dia- 
gram depicted in Ref . , we see that the softcore phase 
diagram is much more fruitful. And on the other hand, 
it is more practical to be verified by experiments because 
U is always finite. Comparing with the softcore bosons 
in a square lattice |l5| , the frustration of a triangular lat- 
tice opens a wider window to explore the SS state. The 
CDW III solid state at integer filling as the counterpart of 
the MI state is another new feature in the present model 
and is worth to be experimentally measured. Mermin- 



Wagncr theorem obstructs the various phases described 
in the present work to be detected in any finite temper- 
ature. However, it may become possible if we trap the 
atoms by an additional shallow trapping potential as the 
Bose-Eistein condensate is observed in a trapped two- 
dimensional cold Bose atom gas. 

This work was supported in part by Chinese National 
Natural Science Foundation. 



[3] 
M 

[5] 



[6 

[7; 

[§ 
[9 
[10 

[11 

[12 
[13 
[14 
[15 
[16 
[17 

[18 

[19 
[20 



N. Penrose and L. Onsager, Phys. Rev. 104,576 (1956) 
A. F. Andreev and I. M. Lifshitz, Sov. Phys. JETP 29, 
1107 (1969); G. Chester, Phys. Rev. A 2, 256 (1970); A. 
J. Leggett, Phsy. Rev. Lett. 25, 1543 (1970). 
E. Kim and M. H. W. Chan, Nature 427, 225 (2004); 
Science 305, 1941 (2004). 

D. M. Ceperley, B. Bernu, Phys. Rev. Lett. 93, 155303 

(2004) ; N. Prokofev, B. Svistunov, ibid. 94, 155302 

(2005) ; D. E. Galli, M. Rossi, L. Reatto, Phys. Rev. B 
71, 140506(R) (2005). 

M. Greiner, I. Bloch, O. Mandel, T. W. Hansch, and 

T. Esslinger , Phys. Rev. Lett. 87, 160405 (2001). C. 

Orzel, A. K. Tuchman, M. L. Fenslau, M. Yasuda, and 

M. A. Kasevich, Science 291, 2386 (2001). M. Greiner, O. 

Mandel, T. Esslinger, T. W. Hansch and I. Bloch, Nature 

415, 39 (2002). M. Greiner, O. Mandel, T. Esslinger, T. 

W. Hansch and I. Bloch, Nature 419, 51 (2002). 

A. Griesmaier, J. Werner, S. Hensler, J. Stuhler, and T. 

Pfau, Phys. Rev. Lett. 94, 160401(2005). 

D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner and P. 

Zoller, Phys. Rev. Lett. 81, 3108 (1998). 

G. G. Batouni and R. T. Scalettar, Phys. Rev. Lett. 84, 

1599 (2000); 

S. Wessel and M. Troyer, Phys. Rev. Lett. 95, 127205 
(2005). 

D. Heidarian and K. Damle, Phys. Rev. Lett. 95, 127206 
(2005). 

R. G. Melko, A. Paramekanti, A. A. Burkov, A. Vish- 
wanath, D. N. Sheng, and L. Balents, Phys. Rev. Lett. 
95, 127207 (2005). 

M. Boninsegni and N. Prokofev, Phys. Rev. Lett. 95, 
237204 (2005). 

S. V. Isakov, S. Wessel, R. G. Melko, K. Sengupta, and 

Yong Baek Kim ,E-preprint cond-mat/0602430 

R. G. Melko, A. Del Maestro, and A. A. Burkov, 

cond-mat/0607501 

P. Sengupta, L. P. Pryadko, F. Alet, M. Troyer, and G. 
Schmid, Phys. Rev. Lett. 94, 207202 (2005). 
G. G. Batrouni, F. Hebert, and R. T. Scalettar, E- 
preprint, |cond-mat/06055i8| 

A. W. Sandvik and J. Kurkijarvi, Phys. Rev. B 43, 5950 
(1991); A. W. Sandvik, J. Phys. A 25, 3667 (1992); A. 
W. Sandvik, R. R. P. Singh, and D. K. Campbell, Phys. 
Rev. B 56, 14510 (1997). 

A. W. Sandvik, Phys. Rev. B 59 R14157701 (1999); O. f. 
Syljuasen and A. W. Sandvik, Phys. Rev. E 66, 046701 

(2002) . 

N. G. Zhang and C. L. Henley, Phys. Rev. B 68, 14506 

(2003) . 

E. L. Pollock and D. M. Ceperley, Phys. Rev. B 36, 8348 
(1987). 



